The evolving distribution of hard partons traversing a hot strongly interacting plasma 
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We present results for the evolution of the momentum distribution of a hard parton traversing a 
brick of quark-gluon plasma, considering both bremsstrahlung and collisional energy loss. The com- 
plete leading order transition rates are included. We find a significant modification of the evolving 
momentum distribution compared to results obtained by employing an approximate implementation 
of the collisional energy loss that uses the mean energy loss, together with momentum diffusion. 



INTRODUCTION 

High transverse momentum jets emerging from the 
central rapidity region in heavy ion collisions can provide 
important information on the created hot quark-gluon 
plasma (QGP). After the discovery of jet quenching at 
the Relativistic Heavy Ion Collider [l|, E| a lot of progress 
has been made toward using jets as a quantitative to 
mographic probe of the QGP 



lets as a quantitative 
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12 . Ila . llJ, lla, llfl, lljj . Jet quenching refers to the sup 
pression of high transverse momentum hadrons, such as 
7r° and rj mesons in central Au + Au collisions compared 
to expectations from measurements in p + p collisions. 
This suppression has been attributed to the energy loss 
of hard px partons due to induced gluon bremsstrahlung 
in the hot quark-gluon plasma phase. Several theoretical 
formalisms have been established to describe the energy 
loss due to bremsstrahlung @, E, El, Efll, El El | • 

In addition, collisional energy loss should play a role 
in the suppression of high momentum jets. However, its 
importance compared to radiative losses has been un- 
der discussion. On the one hand, there have been ar- 
guments for that its contribution is small compared to 
radiative energy loss. For example, early estimates us- 
ing asymptotic arguments indicated that the radiative 
energy loss is much larger than the elastic energy loss 
[23l ]. Furthermore, in 21 1 radiative energy loss in the 
light-cone path integral approach and collisional energy 
loss employing the Bjorkcn method were compared and 
collisional energy loss was found to be smaller. In 24 1 
phenomenological limits on radiative vs. collisional en- 
ergy loss were derived by considering quadratic vs. linear 
pathlcngth dependence, again finding that any elastic en- 
er gy lo ss component has to be small. On the other hand, 
in |25l. Efl it was found that collisional energy loss has a 
significant influence on jet quenching, and recent studies 



27, [28| also point in this direction. See in addition Refs. 
2M and [aT 



In several works 13, HI the evolution of the jet mo- 
mentum distribution including both radiative and colli- 
sional energy loss has been computed. In both works it 
was found that, although the total energy loss of light 
quarks is dominated by radiative processes, the shape 
of the momentum distribution depends significantly on 



whether elastic processes are included or not. However, 
motivated by the fact that collisional energy loss is dom- 
inated by small momentum transfers, in both [l2j and 
[IH elastic collisions were approximated by a mean en- 
ergy loss, accompanied by momentum diffusion. 

Because it was found that the evolving shape of the 
jet momentum distribution can be strongly affected by 
including elastic processes, it is important to go beyond 
this approximation, which from here on we refer to as 
diffusion method. In [28| modifications to this method 
due to fluctuations have been investigated and found to 
be important. In this work we include the full pertur- 
bative transition rates to obtain a more precise result 
and test the validity of the approximation. We simulate 
the situation of a high momentum parton traversing a 
brick of static QGP, considering collisional and radiative 
energy loss 0, EM EH E2 separately and combined. In 
both cases we study the evolution of the parton's entire 
momentum distribution, comparing the diffusion method 
to the direct leading order calculation. 



FORMALISM 

The jet momentum distribution P(E,t) = 
dN(E,t)/dE evolves in the medium according to a 
set of coupled Fokker-Planck type rate equations of the 
form Q: 



dP(E) 
dt 



= I dw (p(E 



+(j) dT(E+u J ,u J )_ p{E) dT(E,u J : 



(i) 



where dT(E,Lu)/dio is the transition rate for processes 
where partons of energy E lose energy uj. The oj < part 
of the integration incorporates processes which increase a 
particle's energy. The radiative part of the transition rate 



is discussed in 0, 3, 33| and the elastic part is discussed 
in detail in the following. 

First, we briefly review the transition rates of the dif- 
fusion method used in 151 ] . It was argued in this work 



that it should be an adequate procedure to approximate 
transition rates due to elastic collisions of the hard par- 
ton with the thermal medium using the mean energy loss 
(drag term) and a momentum diffusion term as dictated 
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by detailed balance. The motivation for this was that 
compared to radiative energy loss, collisional losses are 
more dominated by small energy transfers because the 
contribution to the mean energy loss rate dE/dt from 
elastic collisions is only logarithmically sensitive to large 
energy transfers. On the other hand, the radiative con- 
tribution is a power-law, for large values of the radiated 
energy. It was thus anticipated that, as long as radiative 
energy losses dominate jet quenching, this method will 
adequately address the effects of elastic collisions. Wc 
investigate this assertion quantitatively in this work by 
going beyond the diffusion approximation and monitor- 
ing the entire profile of the energy spectrum. 

In practice, the analytic results for the energy loss rate 
dE/dt HI 0, H3 were incorporated into |T) by in- 
troducing the drag term, (dE / dt)dP(E) / dE , and the dif- 
fusion term, T{dE/dt)d 2 P(E)/dE 2 . Then Eq. ^ was 
discretized, such that J dtu — ► Acj^ w=nAw , and 

T(E + Alu,Alo) = (1 + / B (Aw))(Aw)- 1 dE/dt, 
T(E,-Au) = /b(Acj)(Acj)" 1 dE/dt, (2) 

which yields the right energy loss rate and preserves de- 
tailed balance, for small enough Auj. Using this method 
the momentum distribution evolves like a Gaussian with 
increasing width, shifted by the mean energy loss at the 
regarded time. 

Next, two methods for computing the transition rates 
beyond the diffusion approximation are discussed. One 
may start from the expression for the transition rate 

^(P T\-rl f f 

du; 1 ' ' } k J (2tt) 3 J (2tt) 3 lGpp'kk' 
x 8(p — p — uj)S(k — k — lu) 
x\M\ 2 f(k,T)(l±f(k',T)), (3) 

where p = E = |p| and p' — |p'| are the absolute val- 
ues of the three-momenta of the incoming and outgoing 
hard parton, respectively, and k = |k| and k' = |k'| are 
those of the incoming and outgoing thermal parton. In 
addition, lu = p — p' = k' — k is the transferred energy. 
The distribution functions / are either Fermi-Dirac or 
Bose-Einstein distributions depending on the nature of 
the thermal parton involved. The + or — sign appears 
accordingly, with — for Pauli blocking and + for Bose en- 
hancement, and dk describes the degeneracy of the ther- 
mal parton. 

One way to calculate the transition rate from Eq. (|3|) 
is to follow the energy loss calculation of Braaten and 
Thoma 34, 3f| by shifting the k' integration to the ex- 
changed momentum q = p p' = k' k, and sepa- 
rately computing the contributions from soft ~ gT and 
hard ~ s/ET momentum exchange. For that matter, in 
[Hi], [3H an intermediate separation momentum q* was 
introduced, which dropped out in the final result, when 
adding the hard and soft contribution of the energy loss 
rate. 



Unfortunately, in the calculation of the transition 
rates, q* will not drop out in the final result, because 
some necessary approximations and simplifications due 
to symmetries of the integrand that could be done in 
the calculation of the energy loss rate cannot be done in 
this case. For example, to maintain detailed balance, the 
exact expression f(k)(l—f(k')) must be kept and, in par- 
ticular, cannot be replaced by (f(k) — f(k'))/2 as done in 
34j. In fact, this procedure would lead to negative tran- 
sition rates. Keeping the full expression /(fc)(l — f{k')) 
in the soft part of the calculation will not lead to a loga- 
rithmic divergence, and the above mentioned cancellation 
between the infrared divergence of the hard part and the 
ultraviolet divergence of the soft part will not occur. Wc 
will show below how this problem can be circumvented 
using the method of calculating collisional energy loss 
presented in [38j | . 

For now, we separately compute the hard and the 
soft contribution introducing the separation scale q* and 
match both parts at an intermediate value of q* 2 ~ 



\fETm g (to be precise we use the g*-value for which 
the analytic results for the soft and hard part of dE/dt 
H, H3, [Si, [13l become equal). m 2 g = 1/2(1 + N f /6)g 2 T 2 
is the thermal gluon mass, the characteristic soft scale of 
the system. We vary q* around the mean value to get 
a feeling for the uncertainty that enters the result. We 
refer to this method as 'method A' in the following. 

Eq. ([3]) can be rewritten by introducing the discussed 
shift of the integration to q and also choosing q to lie on 
the z-axis, while placing p in the ir-z-planc. Taking the 
limit p — > 00 leads to 



dT 

duj 



(E, lo, T) 



dq 



dk9{q-\to\) 



2k 



2tt 



(2ir)* 16 E 2 J 

-\M\ 2 f(k,T)(l±f(k',T)), 



(4) 



where <t> pq \k q is the angle between the pxq and the k x q 
plane. The integration limits and the (9-function take care 
of the kinematic restrictions q < k + k' , q < p + p' and 
— q < lu < q. 

For E 3> T the scattering amplitude is dominated by t- 
channcl exchange processes, for which the squared matrix 
elements read 



\M\\ q 



9 9 



s 2 + u 2 



\M\ 2 



\M\ 



su 



\M\ 



19 



2g 4 

9 , 



su\ 



17 



su 



, (5) 



with the Mandelstam variables 

S = -^{[(p + p')(k + k') + q 2 } 

q \ kq W{4pp' + t){4kk' + t)} , 



j. 2 2 

t = U! -q . 

u = — s — t . 



(G) 
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k! 



which in the limit of small u> and q becomes 



FIG. 1: Feynman diagram for the amplitude that contributes 
to collisional energy loss in a QCD medium. The hard thermal 
loop resummed gluon propagator is indicated by a gray blob. 



I-Mlfj, where i £ {q, g} is the squared matrix element 
for the process of a hard parton of type i scattering with 
a thermal parton of type j. Plugging these matrix ele- 
ments into Eq. ([4]) and solving the k- and ^-integrations 
numerically after imposing q* as a lower cutoff on the 
g-intcgration, we obtain the hard part of the transition 
rate. 

To obtain the soft part, we can calculate the matrix 
element for e.g. qq — > qq scattering from the diagram 
shown in Fig.[T]using the effective thermal gluon propaga- 
tor. This method was used in 3J], and further validation 
was put forward in [38| . The matrix element reads 



\M\ 2 qq =-<? 4 Tr(f VW^WV) 

xD^(Q)Dt p (Q), (7) 



where four-momenta are denoted by capital letters, e.g. 
Q = (w,q). The effective thermal gluon propagator in 
the Coulomb gauge is given by 

D^(Q) = S^°6 u0 A L (uj 7 q) + P£ v A T (u,q) , (8) 

where PjP = 0, Pj? = <5' J — <f g- 7 is the transverse projec- 
tor, and 



A £ (w,g) 
At(w, q) 



-1 



q- - mj. 



zln(§±±) - 2 



-1 



q 2 (x 2 — 1) — m 2 



|(1 _ a?) In (=±1) 



are the longitudinal and transverse gluon propagators, 
with x =bj/q. We find 



2tt 

8 - 9 y{[(k + k>) 2 - q 2 ]\A L \* 



1 / 



[(fc + fc') 2 + 9 2 ]|AT| 2 }, (9) 



2- 



(10) 



Note that we keep a factor of k' = k + u) even in this 
limit to fulfill detailed balance exactly. Inserting this ex- 
pression into Eq. ((4| and numerically solving the integrals 
after imposing q* as an upper cutoff on the g-integration, 
we obtain the final result for the soft part. Results for the 
other processes follow analogously and only differ by the 
appropriate prcfactors. Transition rates for gluon energy 
loss are obtained by multiplying those for quarks by 9/4 
[3^ . The complete result is obtained by adding the hard 
and soft contributions. 

Alternatively, we can use the same method as used for 
the soft part above and apply it to the whole transferred 
momentum range as done in [3Sl |. In this case we do 
not use approximations assuming small to or q as done 
in [Hi and our Eq. (flQ|). but insert Eq. © directly into 
Eq. ([J| . Again, we solve the integrals in Eq. numeri- 
cally and now obtain ^-independent transition rates that 
interpolate smoothly between the soft and hard regime of 
method A. We will refer to this method as 'method B' in 
the following. Note that because we are interested in the 
transition rate dT/du {E,ui,T) and not only in the in- 
tegrated energy loss dE/dt (E,T), in Eq. Q we can not 
replace the Bose enhancement and Pauli blocking terms 
(1±/0',T)) by 1 as done in 

In Fig. [2] we compare the transition rates computed 
using the two different methods described, and find a very 
good agreement between the two results as long as we 
choose q* to lie at the geometric mean between the hard 
~ V ET and the soft ~ gT scale in method A. The band 
in Fig.[2]indicates the variation of the result when varying 
q* around that mean value by a factor of two in both 
directions. The slight difference for larger u> stems from 
the fact that in method B it was implicitly assumed that 
—t<^s and hence s 



34j and [36l | ) , which was 



-u (see 

not done in our calculation of the hard part in method A 
(see Eq. ([5])). We have verified that when using s = — u in 
the calculation of the hard part in method A, the results 
from method A agree with those from method B at large 
lo. The difference at large ui will be negligible for the jet 
evolution because elastic energy loss is dominated by soft 
momentum transfers. 

As in [HI, we include the conversion processes (quark 
- anti-quark annihilation, pair production, and Compton 
scattering) that turn a quark into a gluon and vice versa. 
There is no logarithmic enhancement of the u-channel 
exchange process as found for tagged heavy quark en- 



ergy loss in [39( and [40(. See also [4l|. Also we neglect 
sublcading constants arising from the u- and s-channel 
exchange. 
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FIG. 2: (Color online) Transition rate at E = 10 GeV and T = 
200 MeV for the process qq — ► qq calculated in two different 
ways. Nf = 3. 



and soft part is always positive. So the sum of both is 
generally larger than the sum of the analytic solutions. 
Hence the difference in the mean energy loss is related 
to the uncertainty due to using perturbativc methods at 
a s = 0.3. Different approximations, all valid at infinitely 
small coupling, and hence all leading to the same result in 
this regime, can lead to different results when extrapolat- 
ing to larger couplings. Because in the diffusion method 
we use the analytic expressions for the energy loss rate 
but do not in the other methods, we find a qualitatively 
similar result to that of the direct energy loss rate calcu- 
lation described above. 

At T = 200 MeV, the difference in the mean energy 
loss after 5 fm is approximately 30% for purely collisional 
energy loss. It will be significantly smaller in the com- 
bined radiative and collisional calculation, because of the 
dominance of the radiative part. 



RESULTS 

We present results for the evolution of the momen- 
tum distribution of a parton with E = 10 GeV, passing 
through "a brick" of QGP at constant temperatures T = 
200 MeV and T = 300 MeV, using a s = g 2 /{Air) = 0.3 
and Nf — 3. The values of these parameters were chosen 
as typical of the range in many phenomcnological appli- 
cations. Results for a parton with E = 100 GeV were 
found to be qualitatively very similar. We will show the 
distributions after the jet has passed through a 2 fm and 
5 fm thick medium. 

First we compare results from the diffusion method 
used in [l5[ to those obtained using methods A and B 
explained above for purely elastic energy loss. Fig. [3] 
shows the momentum distributions after the quark has 
passed 2 fm of a T = 200 MeV plasma, Fig.H those af- 
ter 5 fm. The results obtained using methods A and B 
agree very well with each other but differ significantly 
from the Gaussian shape emerging using the diffusion 
method. The same qualitative result is obtained for a 
higher temperature T = 300 MeV. 

The total energy loss for both methods A and B is 
slightly larger than for the diffusion method. To under- 
stand this, let us consider the direct calculation of the 
energy loss rate. In 42], 43^ . where the energy loss rate 
of a heavy fermion was calculated numerically and com- 
pared to the analytic result by Braaten and Thoma, it 
was found that the numerically calculated energy loss 
rate was usually larger than the analytic solution. Only 
for very small coupling the two results agreed around an 
intermediate q* . The reason for the difference is that the 
analytic solution involves further approximations that al- 
low both the hard and the soft contribution to become 
negative (which is necessary for the cancellation of q* in 
their sum), while the numerical result for both the hard 
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FIG. 3: (Color online) Jet momentum distribution at T = 
200 MeV for an initial quark jet with E = 10 GeV after 2 fm, 
for a 3 — 0.3. The comparison shows the result for collisional 
energy loss obtained by the diffusion method and method A 
and B described above. 



Next, we show results including radiative energy loss. 
To get a feeling for how much the variation of q* in 
method A affects the final result, in Fig. [5] we show the 
jet momentum distribution as found after the parton has 
passed a static T = 300 MeV QGP of length 2 fm, using 
q* at the discussed mean value between the hard and soft 
momentum scale and for twice and half that value. The 
band shown is a measure of the uncertainty related to 
the (^-dependence in method A. 

Figs. [6] to [9] show the results including both radiative 
and collisional energy loss for two different temperatures 
and two different path lengths through the medium. The 
result for purely radiative energy loss is shown as well. 
It can be seen that also in the combined radiative and 
collisional calculation the differences between the results 
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FIG. 4: (Color online) Jet momentum distribution at T = 
200 MeV for an initial quark jet with E = 10 GeV after 5 fm, 
for a s = 0.3. 
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FIG. 5: (Color online) (/"-dependence in method A and com- 
parison to method B and the diffusion method. The result of 
method A using the geometric mean between the hard and 
the soft scale for q* coincides with the result obtained using 
method B. 



FIG. 6: (Color online) Jet momentum distribution at T — 
200 MeV for an initial quark jet with E = 10 GeV after 2 fm, 
for a s = 0.3. The comparison shows the result for radiative 
energy loss only, collisional and radiative energy loss using 
the diffusion method, and the ones obtained by using method 
A and B described above. 




FIG. 7: (Color online) Jet momentum distribution at T = 
200 MeV for an initial quark jet with E = 10 GeV after 5 fm, 
for a s — 0.3. 



obtained using the diffusion method and those obtained 
using method A or B are significant. 

CONCLUSION 

We studied the evolution of the momentum distribu- 
tion of a high momentum parton traversing a brick of 
static QGP, considering radiative and collisional energy 
loss. For the collisional energy loss we used the lead- 
ing order perturbative expression, which is an improve- 



ment over the method of approximating it by a mean 
energy loss accompanied by momentum diffusion, which 
has been used in recent works [IH, EH ■ 

While the mean energy loss is only slightly larger when 
using the more precise transition rates as opposed to 
the diffusion method, a significant difference in the final 
shape of the momentum distribution is found. It is thus 
entirely conceivable that, depending on the specific ob- 
scrvables, this shape difference could play a large role in 
the interpretation of experimental measurements. Work 
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FIG. 8: (Color online) Jet momentum distribution at T = 
300 MeV for an initial quark jet with E = 10 GeV after 2 fm, 
for a s = 0.3. 
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FIG. 9: (Color online) Jet momentum distribution at T = 
300 MeV for an initial quark jet with E = 10 GeV after 5 fm, 
for a s = 0.3. 



in this direction is ongoing. 

We therefore conclude that it is important to go 
beyond the diffusion approximation, particularly when 
comparing different theoretical energy loss formalisms in- 
cluding both radiative and collisional energy loss. 
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